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1 Abstract 

In this article we report the release of a new program for calculating the emissivity 
of atomic transitions. The program, which can be obtained with its documentation 
from our website www.scienceware.net, passed various rigorous tests and was 
used by the author to generate theoretical data and analyze observational data. 
It is particularly useful for investigating atomic transition lines in astronomical 
context as the program is capable of generating a huge amount of theoretical data 
and comparing it to observational list of lines. A number of atomic transition 
algorithms and analytical techniques are implemented within the program and can 
be very useful in various situations. The program can be described as fast and 
efficient. Moreover, it requires modest computational resources. 


2 PROGRAM SUMMARY 


4 


2 Program Summary 

Title of the program: Emissivity 
Type of program: command line 
Programming language: C++ 

Number of lines: 3110 

Computer: PC running Linux or Windows 

Installation: desktop 

Memory required to execute: case dependent 
Speed of execution: case dependent (~ few seconds) 

Has code been vectorized?: Yes 
Has code been parallelized?: No 
Compilers tested: g++, Dev-C++ 

Number of warnings: 0 

Programming methodology: procedural with object oriented 
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3 Theoretical Background 

In a thermodynamic equilibrium situation an excited atomic state is populated by 
recombination and radiative cascade from higher states, and depopulated by au- 
toionization and radiative decay to lower states. Many recombination lines arise 
from radiative decay and subsequent cascade of strongly autoionizing resonance 
states near the ionization limit. These lines are dominated by low temperature di- 
electronic recombination. The populations of such resonance states are determined 
by the balance between autoionization and radiative decay. When autoionization 
dominates, the populations are then given by the Saha equation for thermodynamic 
equilibrium 


where A E t is the energy of the recombined electron in the X^ n _i) + state relative to 
the ionization threshold, and the other symbols have their usual meaning as given 
in Nomenclature § 8. 

The Saha equation, which describes the ratio of different stages of ionization, 
is based on the assumption of Local Thermodynamic Equilibrium (LTE) in a gas 
where collision dominates other physical processes. Consequently, the local velocity 
and energy distributions of particles are given by the Maxwell and Boltzmann 
distributions respectively and a temperature can be defined locally. The Saha 
equation is therefore strictly applicable only if elastic collisions are responsible 
for establishing the energetic distribution of particles. In many practical cases, 
however, atomic processes such as radiative transitions or dynamic effects are more 
important than elastic collisions and the assumption of LTE is not justified within 
the whole energy range. In these cases explicit detailed equilibrium calculations 
are required to determine the velocity and energy distributions of particles over the 
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various energy levels [1], 

To measure the departure of the state from thermodynamic equilibrium a de- 
parture coefficient, b u , is defined as the ratio of autoionization probability to the 
sum of radiative and autoionization probabilities. The value of this coefficient is 
between 0 for radiative domination and 1 for thermodynamic equilibrium. A large 
value of b u (~ 1) is then required to justify the assumption of a thermodynamic 
equilibrium and apply the relevant physics. 

In thermodynamic equilibrium (TE) the rate of radiationless capture equals the 
rate of autoionization, giving [5] 


NMe, = n°ti ( 2 ) 

where N e and iVj are the number density of electrons and ions respectively, g c is 
the recombination coefficient for the capture process, N‘J is the Saha population of 
the doubly-excited state and T“ is the autoionization transition probability of that 
state. In non-TE situation, the balance is given by 

NeNte, = JV„(U + U) (3) 


where N u is the non-TE population of the doubly-excited state and T r u is the 
radiative transition probability of that state. The departure coefficient b u is a 
measure of the departure from TE, and hence is the ratio of the non-TE population 
to the Saha population. Comparing Equations (2) and (3) gives 
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4 Long Write-up 

‘Emissivity’ code is a command line program that was developed to implement the 
atomic emissivity model. Its main functionality is to calculate the emissivity of the 
transition lines from electronic recombination and cascade decay all the way down 
to the ground or a metastable state, that is all free-free, bound-free and bound- 
bound transitions. ‘Emissivity’ is written in C++ computer language and mixes 
procedural with object oriented programming. It consists of about 3000 lines of 
code. The program was compiled successfully with no errors or warnings using 
Dev-C++ compiler on Windows, and g++ compiler on Cygwin and several Linux 
distributions such as Red Hat Enterprise. Sample results produced by these three 
versions were compared and found to be identical. Elaborate internal checks were 
carried out in all stages of writing and debugging the program and the output 
was verified. Thorough independent checks on sample emissivity data produced by 
‘Emissivity’ were performed and found to be consistent. 

The program reads from plain text input hies and writes the results to a main 
plain text output hie. Other secondary output hies can also be produced for par- 
ticular purposes when required. The necessary input hies are 

• The main hie to pass the parameters and inform the program of the required 
calculations. The main parameters in this hie will be outlined in § 5.1. 

• A hie for passing the free-states data and oscillator strengths (/-values) for 
free-free and bound-free transitions with photon energies. The free-states 
data include an index identifying the resonance, its energy position, width, 
configuration, term, 2 J, parity, and a hag marking the energy position data 
as experimental or theoretical. 

• A hie containing the bound-states data. These data include an index iden- 
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tifying the state, its energy, effective quantum number, configuration, term, 
and a flag marking the energy data as experimental or theoretical. 

• A hie containing the oscillator strengths (/-values) for the bound-bound tran- 
sitions. 

The last hie is imported from the R-matrix code [2, 4] output while the rest are 
user made. Two other input data hies are also required if comparison and analysis of 
observational data are needed. One of these is a data hie that contains astronomical 
observations while the other includes mapping information of observational lines 
to their theoretical counterparts. The astronomical data hie contains data such 
as observed wavelength, observed hux before and after correcting for reddening 
and dust extinction, ionic identification, laboratory wavelength, multiplet number, 
lower and upper spectral terms of transition, and statistical weights of lower and 
upper levels. The mapping hie contains the indices of the observational lines and 
their theoretical counterparts. 

In this section we summarize the theoretical background for emissivity calcula- 
tions as implemented in the ‘Emissivity’ code 

• The program starts by obtaining the radiative transition probability T r ul for 
all free-free transitions as given by 

_ a 3 E 2 p gJ u i 

ul igwo ( } 

where the photon energy E p is in Ryd, and the other symbols are defined 
in Nomenclature § 8. This is followed by obtaining the radiative transition 
probability T r ul for all free-bound transitions, as in the case of free-free tran- 


sitions. 
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• The total radiative transition probability T^ for all resonances is then found. 
This is the probability of radiative decay from an upper resonance state to 
all accessible lower resonances and bound states. This probability is found 
by summing up the individual probabilities T r ul over all lower free and bound 
states / for which a transition is possible according to the electric dipole rules; 
that is 

n = Y. c ( 6 ) 

i 

• The departure coefficient, b u , for all resonances is then obtained 


b 


U 


T a 

U 


+ T r 

1 u 


( 7 ) 


where T“ and T^ are the autoionization and radiative transition probabilities 
of state u, and T a is given by 
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• The next step is to calculate the population of resonances by summing up 
two components: the Saha capture, and the radiative decay from all upper 
levels. In thermodynamic equilibrium (TE) the rate of radiationless capture 
equals the rate of autoionization, giving 

KN<e. = N?Ti (9) 


where N e and iVj are the number density of electrons and ions respectively, 
q c is the recombination coefficient for the capture process, N t s is the Saha 
population of the doubly-excited state and T“ is the autoionization transition 
probability of that state. In non-TE situation, the population and depopula- 
tion of the autoionizing state due to radiative decay from upper states and to 
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lower states respectively should be included, and hence the balance is given 
by 

JV«JV ie „ + £jVJ^ = W,(r? + r7) ( 10 ) 

u 

where Ni is the non-TE population of the doubly-excited state, T[ is the 
radiative transition probability of that state, T r ul is the radiative transition 
probability from an upper state u to the autoionzing state Z, and the sum is 
over all upper states that can decay to the autoionzing state. Combining (9) 
and (10) yields 

Jvfr? + £ Jv„ U, = JV,(r? + r[) (ii) 

u 

On manipulating (11) the following relation can be obtained 
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where bi is the departure coefficient of the autoionizing state. This last rela- 
tion is used in calculating the population. 


• The next step is to calculate T() z for the bound-bound transitions. In these 
calculations the /-values can be in length form or velocity form, though the 
length values are usually more reliable, and hence ‘Emissivity’ reads these 
values from the /-values hie produced by the R-matrix code. This is followed 
by Ending T„ for the bound states by summing up T r ul over all lower bound 
states /, as given earlier by (6) for the case of resonances. 

• The population of the bound states is then obtained 
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where u includes all upper free and bound states. 

• Finally, all possible free-free, free-bound and bound-bound transitions are 
found. The emissivity of all recombination lines that arise from a transition 
from an upper state u to a lower state l is then computed using the relation 

Sul = N u V r ul hu (14) 

where v is the frequency of the transition line. 

• Apart from the normal debugging and testing of the program components 
to verify that they do what they are supposed to do, two physical tests are 
incorporated within the program to validate its functionality and confirm 
that no serious errors have occurred in processing and producing the data. 
These tests are the population-balance test and the metastablc test. The 
first test relies on the fact that the population of each state should equal the 
depopulation. This balance is given by the relation 


E N i rr p = v E r « < 15 ) 

j>i k<i 

The second test is based on the fact that the total number of the electrons 
leaving the resonances in radiative decay should equal the total number ar- 
riving at the metastable states. This balance is given by the relation 

E N i r i = E ^ < 16 ) 

Vj Vi,k>i 


where i is an index for metastable states and j is an index for resonances. 
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5 Input and Output 

The detailed technical description of the program input and output files is given 
in the program documentation which can be obtained from www . scienceware . net 
website. However, in this section we give a general description of the nature of the 
input parameters and the expected output data so that the user can decide if the 
program is relevant for the required purpose. 

5.1 Input 

Various physical and computational parameters are required for successful run. 
Some of these are optional and are only required for obtaining particular data. 
The main input parameters are 

• The temperature in K and the number densities of electrons and ions in 
m -3 . The program uses a temperature vector to generate data for various 
temperatures. These values can be regularly or irregularly spaced. 

• The residual charge of the ion, which is required for scaling some of the data 
obtained from the R-matrix code output hies. 

• Boolean hag for reading and processing observational astronomical data ob- 
tained from an external input data hie. A mapping hie is also required for 
this process. If the user chooses to read astronomical data, the data lines will 
be inserted between the theoretical lines in the main output hie according to 
their observed wavelength. 

• Boolean hag to run the aforementioned physical tests. 

• Parameters for generating and writing normalization emissivity data for the- 
oretical and observational lines to the main output hie alongside the original 


5.1 Input 


13 


emissivity data. The normalization can be internal using the emissivity of one 
of the transition lines produced by the program. It can also be with respect to 
an outside set of emissivity values corresponding to the input temperatures. 
Another possibility is to normalize with respect to an outside set of emissiv- 
ity data in the form of effective recombination coefficients corresponding to 
the specified input temperatures. In this case the wavelength of the line to 
be normalized to is required. The effective recombination coefficient ^(A) is 
defined such that the emissivity e(A) of a transition line with wavelength A 
is 

h c 

e(A) = N e NiQ f (X)— (J.m -3 .s -1 ) (17) 

Another choice for normalization is to be with respect to an outside set of 
emissivity data in the form of effective recombination coefficients correspond- 
ing to a set of temperatures that may be different from those used in the 
input. In this case, the recombination coefficients corresponding to the input 
temperature values are obtained by interpolation or extrapolation using a 
polynomial interpolation routine. 

• Parameters to control the production and writing of the effective recombina- 
tion coefficients, g f , which are equivalent to the given cmissivities. Vacuum 
wavelength of the transition lines will be used due to the restriction on the 
air wavelength formula (A > 2000A). 

• Parameters to control the algorithm to minimize the sum of least square 
deviations of emissivity between the observational lines and their theoretical 
counterparts over the input temperature range. This sum is given by 

5= V ( 4 - 4) 2 

V lines 


( 18 ) 
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where e ° and e f are the normalized observational and theoretical emissivi- 
ties respectively. This algorithm also offers the possibility of computing the 
temperature confidence interval and its relevant parameters. The individ- 
ual square differences for each one of the least squares and their percentage 
difference which is given by 


( e N 0 - ef 

\ $ 


x 100 


(19) 


can also be written to the main output hie. The confidence interval for the 
goodness-of-ht index y 2 can also be found if required by providing the degrees 
of freedom rj and the change in the goodness-of-ht index Ay 2 . 

• Parameters to run and control an algorithm for finding the decay routes to 
a particular state, bound or free, from all upper states. As the number of 
decay routes can be very large (millions and even billions) especially for the 
low bound states, a parameter is used to control the maximum number of 
detected routes. The results (number of decay routes and the routes them- 
selves grouped in complete and non-complete) are written to an output hie 
other than the main one. The states of each decay route are identified by 
their configuration, term and 2 J. For resonances, the Saha capture term and 
the radiative decay term are also given when relevant. These temperature- 
dependent data are given for a single temperature, that is the hrst value in 
the input temperature vector. 


5.2 Output 

The program produces a main output hie which contains the essential emissivity 
data. The hie starts with a number of text lines summarizing the input data used 
and the output results, followed by a few commentary lines explaining the symbols 
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and units. This is followed by a number of data lines matching the number of tran- 
sitions. The data for each transition includes an index identifying the transition, 
status (FF, FB or BB transition), the experimental state of the energy data of the 
upper and lower levels, the attributes of the upper and lower levels (configuration, 
term, 2 J and parity), wavelength in vacuum, wavelength in air for A > 2000A, the 
emissivities corresponding to the given temperatures, the normalized emissivities 
and the effective recombination coefficients corresponding to the given emissivities. 
Writing the normalized emissivities and the effective recombination coefficients is 
optional and can be turned off. As mentioned earlier, if the user chooses to read 
the astronomical data, the observational lines will be inserted in between the the- 
oretical lines according to their lab wavelength. 

There are other subsidiary output hies. One of these is a hie containing the 
results of the minimization algorithm. This hie contains information on r/, Ay 2 , 
a, temperature at minimum y 2 , and the temperature limits for the confidence 
interval. This is followed by the temperature array with the corresponding least 
squares residuals, the ay 2 and the y 2 arrays as functions of temperature. Another 
output hie is one containing mapping data used for testing purposes. A third 
output hie is a decay routes hie which was outlined previously. 
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8 Nomenclature and Notation 

a fine structure constant (= e 2 /(hc4ne 0 ) ~ 7.2973525376 x 10~ 3 ) 

T a autoionization transition probability (s _1 ) 

r r radiative transition probability (s^ 1 ) 

T r ul radiative transition probability from upper state u to lower state l (s _1 ) 

Ay 2 change in goodness-of-fit index 
A r width of resonance (J) 

A E energy difference (J) 

e Q permittivity of vacuum (~ 8.854187817 x 10 _12 F.m _1 ) 

£ cmissivity of transition line (J.s -1 .m -3 ) 

e u i cmissivity of transition line from state u to state l (J.s^hm -3 ) 
e g normalized observational emissivity 

Et normalized theoretical emissivity 

rj number of degrees of freedom 

A wavelength (m) 

v frequency (s -1 ) 

g c recombination coefficient for capture process (m 3 .s _1 ) 
g f effective recombination coefficient (m 3 .s _1 ) 

t 0 atomic time unit (= H/E h ~ 2.418884326505 x 10~ 17 s) 

X 2 goodness-of-fit index 

a scaling parameter in y 2 conhdence interval procedure 
A angstrom 

b departure coefficient 

b u departure coefficient of upper state 

c speed of light in vacuum (299792458 m.s -1 ) 
e elementary charge (~ 1.602176487 x 10 -19 C) 
e~ electron 

E energy (J) 

Eh Hartree energy (~ 4.35974394 x 10~ 18 J) 

E p photon energy (J) 

/ oscillator strength 

f u i oscillator strength of transition between upper state u and lower state / 
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9 

h 

h 

J 

k 

L 

m 

m e 

N 

N e 

Ni 

N s 

R 

r 

S 

t 

T 

T 

X n -\- 

V( n _i)+ 


statistical weight in coupling schemes (= 2J + 1 for IC) 

Planck’s constant (~ 6.6260693 x 10 -34 J.s) 

reduced Planck’s constant (= h/2-n ~ 1.0545717 x 10 -34 J.s) 

total angular momentum quantum number 

Boltzmann’s constant (~ 1.3806505 x 10 -23 J.K -1 ) 

total orbital angular momentum quantum number 

mass (kg) 

mass of electron (~ 9.10938215 x 10~ 31 kg) 
number density (m -3 ) 
number density of electrons (m -3 ) 
number density of ions (nr -3 ) 

Saha population (m -3 ) 

resonance matrix in R-matrix theory 

radius (m) 

total spin angular momentum quantum number 
time (s) 

temperature (K) 

electron temperature (K) 

ion of effective positive charge n 

bound state of recombined ion with effective positive charge n — 1 


Abbreviations 

BB Bound-Bound 
C Coulomb 
F Faraday 

FB Free-Bound 
FF Free-Free 
IC Intermediate Coupling 
J Joule 

K Kelvin 
kg kilogram 

LTE Local Thermodynamic Equilibrium 
(.)/ lower 
m meter 
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Ryd Rydberg 

s second 

TE Thermodynamic Equilibrium 

(■)« upper 

(.) u i upper to lower 

V for all 

Note: units, when relevant, are given in the SI system. Vectors are marked with 
boldface. Some symbols may rely on the context for unambiguous identification. 
The physical constants are obtained from the National Institute of Standards and 
Technology (NIST) website [3]. It should be remarked that most subscripts (e.g. I 
and u ) are dummy variables, that is they are subject in their interpretation to the 
context and are not specific to the particular cases they are referring to. 


